»*wswwn igwMMfcawwjawMW  minpwammupwi  m^wwwpwimbwi*  w«ww*t»  -■ 


HJCHANIC3  G?  GRANULAR  MEDIA  * 


For  a number  of  years,  there  has  been  under  development  a mtbemtlcal 
theory  of  the  mechanical  behavior  of  materials  composed  of  discrete  elastlo 
grains  In  direct  o on  tact.  Eventually,  the  theory  ie  intended  to  predict 
stress-strain  relational  atreae  distributional  vibrations,  wave  propagation 
phenomena  and  orltexia  of  failure  for  such  aateriala  as  are  found  in  a bed  of 
dry  sand  or  the  pile  of  grains  in  the  carbon  microphone.  The  line  of  attaok, 
which  has  been  the  aoet  fruitful,  begins  with  a consideration  of  the  local 
forcea  and  deformations  at  the  oostaot  surfaces  between  adjacent  grains* 
Because  of  the  extraordinarily  complex  nature  of  the  problem  the  grains  have 
been  idealised  as  like  spheres  in  regular  arrays.  Even  with  this  simplifica- 
tion, at  least  until  recently,  only  the  component  of  fores  normal  to  each 
contact  surface  has  been  taken  into  aocount  [1,2, 3, 4],  The  relations  between 
normal  force,  N , contact  radius,  a , and  displacement,  06  , are  obtained  fren 
the  Hert*  theory  of  contact  of  elastic  bodies  [5]* 
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where  R is  the  radius  of  the  spheres,  V is  Poisson* s ratio  and^M  the  shear 
modulus  of  the  material  of  the  spheres.  Of  special  interest  is  the  normal 
compliance 
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* lecture  presented  at  the  Second  U.S.  National  Congress  of  Applied  Mechanics, 
Ann  Arbor,  Michigan,  on  June  16,  1954. 
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The  non-linearity  of  these  relatione  gives  the  first  inkling  of  dynamical 
difficulties  in  addition  to  the  purely  geometrical  on es.  The  behavior  of  a 
granular  mate rial  may  be  expected  to  depend  strongly  on  the  initial  stress 
which,  in  turn,  affeots  the  role  of  the  elAstio  constants  of  the  individual 
grains.  The  early  forms  of  the  theory  predict  wave  velocities  proportional 
to  tfas  sixth  root  of  an  initial  isotropic  pressure  and  the  cube  root  (rather 
than  the  usual  square  root)  of  the  shear  modulus  of  the  grain  [2,3,4]*  These 
relations  have  been  confirmed  experimentally  [2,13]  but  absolute  velocities 
are  In  poor  agreement  when  the  theory  does  not  include  the  effect  of  tangen- 
tial components  of  force  s tween  grains c Zt  is  the  purpose  of  this  lscture 
to  dleouss  seme  of  the  js  o blame  and  consequence#  of  including  considers tloa 
of  tangential  force*  ».c  the  ecataot  surfaces* 

Corresponding  to  the  Barts  theory,  there  ia  a solution  of  the  equations 
of  elasticity  [6,7]  which  takes  into  account  a mono  tool  cal  ly  Increasing  tangen- 
tial force  ( T ) subsequent  to  the  application  of  a normal  force.  Zt  is  found 
that  a tangential  force,  no  matter  hew  small,  produces  infinite  tangential 
traction  (X  ) at  the  edge  of  the  contact  surface  (see  Fig.  1)  if  it  is  assumed 
that  there  ia  no  relative  displacement  of  opposing  points  on  the  contact  sur- 
face. Accordingly,  it  is  assumed,  in  the  theory,  that  such  a relative  dis- 
placement does  take  place  and,  because  of  symmetry,  it  occurs  on  an  annulus. 
Rirther,  the  outer  edge  of  the  annulus  is  assumed  to  coincide  with  the  edge  of 
the  o on tact  surface,  because  it  is  there  that  the  infinite  traction  would 
otherwise  occur.  The  boundary  conditions  of  the  theory  of  elasticity  require 
that  there  be  specified,  on  the  annulus,  the  tangential  traction  or  displace- 
ment or  a relation  between  the  two.  In  this  case  it  has  been  assumed  that 
the  tangential  component,  t*  , of  traction  at  each  point  of  the  annulus  is 
proportional  to  the  normal  a exponent,  O'  , at  that  point,  Fhysically , this 
is  to  say  that  slip  takes  placa  on  ths  annulus  in  such  a vay  that  Coulomb* j 


lav  of  friction  hold*  at  each  point,  T-fO"  where  is  the  Berts 

noraal  pressure  and  f Is  a constant  coefficient  of  friotion*  The  resulting 
distribution  of  tangential  traction  over  the  entire  contaot  surface  Is  Ulus-* 
t rated  in  fig,  1. 

As  the  tangential  force  is  increased,  the  theor7  predicts  that  the 
inner  radius  ( c ) of  the  annulus  of  slip  diminishes  according  to  the  law 


(4) 


At  the  ea.ee  time,  the-  relative-  tangential  dlsplaoeeent,  5 , of  distant  points 
in  the  two  spheres  depends  on  the  tangential  force  according  to 
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This  relation  is  shown  in  Fig*  2 along  with  experimental  data,  obtained  by 
Johnson  [8]  with  steel  spheres,  which  confirm  it*  The  tangential  coapllanoe, 
to  be  c capered  with  Equation  (j).  It 


The  next  step,  in  the  study  of  local  effects  at  the  contaot  surfaces, 
was  to  determine  the  consequences  of  reversal  of  the  sense  of  the  tangential 
force  [9]*  If  the  tangential  force,  after  reaching  a magnitude  T^T* 
is  diminished,  the  force-displacement  relation  is 
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This  relation  is  shown  as  the  curve  PRS  in  Fig*  3.  Kara  a new  complication 
is  seen  to  enter,  namely,  the  inelastic  (as  distinguished  from  non-linear 
elastic)  character  of  the  tangential  load-displacement  relation.  In  the  case 
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where  the  tangential  force  oscillates  between  & T*  (where  jr*|  < f N ), 
three  Important  conoluaione  were  reached  t (l)  slip  is  confined  to  an  annulus 
whose  Inner  radius  Is  given  by  Equation  (4)  'dth  c and  T replaced  by  c* 
and  T*j  (2)  the  amplitude  of  the  relative  displacement  of  the  spheres  la 
given  by  Equation  (5)  with  T replaced  byT*  j (3)  the  force-displaceaent 
curve  is  a loop  (fig*  3)  enclosing  an  area  which  represents  the  energy  dis- 
sipation per  cycle  1 
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ill  of  these  conclusions  have  been  subjected  to  experimental  test. 

Tests  by  Mindlln , Mason , Osmer  and  Deresiewics  [10]  were  made  with  a 
pile  of  three  polished  glass  lenses , pressed  together  with  a normal  force 
following  which  an  oscillating  transverse  force  was  applied  to  the  central 
lens  at  60  c.p.s.  (Fig.  4)»  According  to  the  theory,  relative  displacement 
at  the  contact  surface  occurs  only  on  an  annulus,  so  that  wear  patterns  should 
be  observed  only  there  and  with  inner  radius  given  by  Equation  (4)*  Such 
patterns  were  observed  (Fig.  5)  and  the  comparison  of  their  dimensions  with 
those  predicted  by  the  theory  Is  shown  in  Fig.  6.  Measurements  were  also 
made  of  energy  dissipation.  At  large  amplitudes  these  conformed  with  Equation 
(8)  but,  at  small  amplitudes,  the  energy  dissipation  varied  as  the  square  of 
the  tangential  force  rather  than  the  cube  as  the  theory  requires.  Biis  watt 
evidence  that  a velocity  dependent  factor  contribute  to  energy  dissipa- 

tion in  addition  to  the  static  considerations  on  which  Equation  (8)  is  based* 
An  extensive  series  of  both  static  and  dynamic  tests  by  Johnson  [8] 
bear  on  many  aspects  of  the  theory*  His  static  experiments  included  leading, 
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unloading,  overloading  and  cyclic  loading*  all  confirming  t ha  behavior  pre- 
dicted b 7 tho  theory.  In  his  dynamic  teste,  conducted  at  4&*5  o.p.s*  with 
a yariaty  of  sphere  diameters  and  normal  leads,  Johnson  obtained  tha  rela- 
tions betvean  tangential  fores  and  displacement  amplitudes  shown  in  Fig*  7* 

As  may  be  seen,  the  theory  is  eery  good  in  this  respect.  Use  sane  series  of 
tests  yltlded  date  on  energy  dissipation  (Tig.  8)  and  in  this  case  the  thmorj 
le  not  satisfactory*  As  nay  be  seen,  in  Fig*  8,  the  energy  dissipation  per 
oyole  at  small  amplitude#  is  again  found  to  vary  as  the  square  of  the  ampli- 
tude, Indicating  the  presence  of  a Telocity  dependent  mechanism  which  cam* 
pletely  overshadows  the  static  mechanisa  at  eery  small  amplitudes*  In  addi- 
tion, there  appears  to  be  e geometrical  faotor,  missing  in  the  theory,  which 
is  important  at  intermediate  amplitudes,  since,  in  that  region , Johnson's 
experiments  reveal  a dependence  of  energy  dissipation  on  both  sphere  diameter 
and  normal  lead,  which  is  not  accounted  for  in  the  theory.  It  is  only  at 
large  amplitudes  (near  gross  sliding)  that  the  theory  appears  to  give  good 

i 

results  for  energy  dissipation  per  oyole*  J 

In  addition  to  normal  and  tangential  forces  on  the  contact  surfaces,  { 

* 

twisting  couples  can  also  be  present  in  a significant  amount  in  cartaln  types 
of  deformation  of  granular  materials.  The  problems  analogous  to  those  described 
above  for  tangential  forces  hare  also  been  solved  for  twisting  couples  [7,11,12]. 

Before  proceeding  to  eseemblagee  of  spheres  it  was  necessary  to  carry 
the  theory  of  pairs  of  spheres  one  step  farther.  Thus  far,  in  both  theory 
and  experiment,  the  normal  force  was  held  constant  during  variation  of  the 
tangential  force.  However,  in  am  assemblage  of  spheres  under  varying  external 
load  or  internal  vibration,  the  normal  and  tangential  forces  on  a single  con- 
tact surface  vary  simultaneously.  In  this  case  the  inelastic  character  of 
the  relation  between  tangential  load  and  displacement  introduces  a very  great 
complication  in  that  it  causes  tha  instantaneous  tangential  forco-displacaaient 


relation  to  depend  on  the  entire  past  history  of  normal  and  tangential  leading* 
Different  phenomena  are  involved  and  different  re stilts  obtained  depending  upon 
vhether  the  nornal  or  the  tangential  force  is  held  constant,  while  the  other 
varleej  whether  they  both  vary,  and  whether  the  sense  of  the  variation  is 
such  that  oas  increases  while  the  other  decreases , both  Increase,  or  both  de- 
crease; whether  their  relative  rate  of  change  is  greater  or  less  than  the 
coefflolent  of  frlotioo;  whether  the  lsnedlate  past  history  of  loading  was 
in  the  same  or  opposite  sense  ss  the  current  loading*  For  example,  suppose 
that,  after  applying  a Dorsal  force  N,  , both  N and  T are  increased  at  an 
arbitrary  relative  rate*  Then,  la  place  of  Equation  (6),  the  tangential,  compli- 
ance is  [9j 
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where  O is  the  instantaneous  radius  of  the  contact  surface.  Compliances  of 
this  type  enter  into  the  prediction  of  failure  loads  of  granular  materials* 

The  implications  of  the  fora  of  Equation  (9)  ara  discussed  below* 

Another  case,  of  interest  in  connection  with  vibrations  of  granular 
naterials,  is  that  in  which,  after  an  initial  normal  force  N«  is  applied, 
the  tangential  force  oscillates  between  drT*  while  the  normal  force  varies 
in  such  a way  that  dti/dT  is  constant*  The  tangential  compliance  during 
the  loading  part  of  the  cycle  is 


where 
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For  the  unloading  port  of  the  070I9  the  signs  of  a and  L are  reversed  in 
Equation  (10).  the  associated  "static"  energy  dissipation  per  oyole  is 


Consider,  now,  a granular  body  coo  posed  of  like  spheres.  If  the  bod/ 

Is  fully  consolidated  the  arrangement  of  the  spheres  is  a faoe-centered  cubic 
or  hexagonal  array,  both  of  these  being  arrangenente  of  densest  packing.  An 
incompletely  consolidated  body  contains  clusters  of  spheres  having  such  pack- 
ing. We  begin  by  considering  an  element  of  a face-centered  cublo  array  of 
spheres  in  equilibrium  under  an  arbitrary  state  of  initial  stress  and  aak 
what  deformation  will  result  frca  an  arbitrary  additional  increment  of  stress. 
This  question  has  been  explored  in  detail  recently  [13] • 

The  elementary  block  of  the  face-centered  cubic  array  is  shown  in 
Fig.  9 , and  the  cceponerts  of  lnoremontal  f orce , d Pij  • , acting  on  it  are 
shown  in  Fig.  1C(#  The  incremental  stress  is  defined  as  the  ratio  of 

the  incremental  forqe  to  the  area  of  a face  of  the  block,  i.e*,  c/0[j  - dP{j/8Rl 
where  R is  the  radius  of  the  spheres.  The  deformation  of  the  block,  resulting 
from  the  application  of  d^ j , can  be  obtained  if  the  increments  of  contact 
force  between  spheres  are  known;  for  then  the  relative  incremental  displace- 
ments of  the  spheres  can  be  found  by  multiplying  by  the  contact  compliances# 
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Each  sphere  In  a face-centered  cubla  array  is  in  contact  with  twelve 

othar  spheres,  Hanca  thara  ora  thirty-six  components  of  contact  forca  on 

aach  sphara,  However,  since  ws  consider,  tempo rari ly , a hoaoganaous  state 

* 

of  Incremental  stress,  eighteen  of  the  components  of  contact  forca  are  equal 
in  pairs,  leaving  only  eighteen  to  be  found,  of  which  six  are  normal  components 
and  twelve  tangential.  Die  latter  are,  in  turn,  related  through  three  equa- 
tions of  moment  equilibrium,  lhe  eighteen  contact  forces  are  related  to  the 
stresses  d<J\j  through  six  independent  equilibrium  equations  so  that,  in  all, 
there  are  only  nine  equations  of  squilifcriua  from  which  to  determine  eighteen 
oootaot  forces;  that  is,  the  problem  la  statically  indeterminate.  It  nay  be 
solved  either  by  introducing  equations  of  compatibility  of  relative  displace- 
ments of  spheres  (there  are  nine  ouch  squat  long)  or  by  starting  with  a set 
of  ouapatible  incremental  strains  d 6y  and  calculating  the  corresponding 
contact  forces.  Die  latter  procedure  is  simpler  since  it  does  not  involve 
the  solution  of  eighteen  simultaneous  equations.  In  either  case  the  incre- 
mental stress-strain  relation  is  found  in  the  form 

3 C‘J>/  d€n  (12) 

Where,  for  the  moat  general  state  of  Jritial  stress,  Cijki  is  • non-symaetrio 
tensor  having  thirty  non-zero  components  when  referred  to  the  principal  axee 
of  the  oubio  array.  These  components  are  linear  functions  of  the  reciprocals 
of  the  eighteen  initial  compliances  associated  with  the  twelve  contact  sur- 
faces, Each  of  the  initial  compliances  depends,  in  turn,  on  the  history  of 
the  initial  stress  according  to  relations  such  as  Equations  (9)  in  which  N 
and  T are  themselves  functions  of  the  stress.  Thus  the  problem  of  solving 
Equation  (12)  to  obtain  a finite  stress-strain  relation  is  a formidable  one 
involving,  as  it  does,  the  solution  of  simultaneous,  ncn-line&r,  integro- 
diffarential  equations.  However,  in  certain  special  cases,  which  can  be 
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rsalijed  in  the  laboratory,  the  integration  of  tba  incremental  stress-strain 
relation  either  can  be  aocoapliahed  or  la  not  necessary. 

An  exaaple  of  a teat  in  which  the  incremental  atreaa-etraln  relatione 
nay  be  uaed  without  Integration  ia  that  of  email  vibration*  In  the  presence 
of  high  initial  a trees.  In  this  case  the  change  In  atrea*  during  vibration 
can  be  made  so  anall  in  c caparison  with  the  initial  stress  that  the  oontaot 
oomplianoee  remain  essentially  constant,  furthermore,  if  the  initial  stress 
is  isotropic  the  incremental  strese-etrmin  relation  reduce*  to  one  of  simple 
oublo  aye  try  with  only  three  ooeffloiantet 
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in  which  (T0  ia  the  initial  isotropic  stress.  In  the  case  of  a high  frequency 
vibration,  c„  f c,*.  and  c**  must  also  have  iaa ginary  parts;  but  the  theory 
is  not  sufficiently  developed  to  write  them  explicitly,  although  Johnson’s 
experiments  give  a good  indication  of  what  their  fora  should  be.  At  present, 
the  imaginary  parts  are  oedtted.  It  is  then  a simple  matter  to  calculate 
wave  velocities  or  frequencies  of  vibration  of  a bar.  Such  bars  were  -construc- 
ted in  tha  following  manner  [13].  A long  rectangular  box,  lined  with  a loose 
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rubber  sheet,  was  carefully  filled  with  3/3"  steal  balls  arranged  in  face- 
centered  cubic  array.  The  *heet  was  thin  folded  over,  sealed  and  evacuated. 
The  external  pressure  locked  the  balls  in  place  so  that  the  solid  "granular 
bar"  could  be  removed  from  the  box  {see  Pig*  11)*  The  balls  were  arranged, 
in  various  bars,  so  that  either  the  [100]  or  the  [110]  direotlon  was  parallel 
to  the  length  of  the  bar  so  as  to  eliminate  ooupling  between  longitudinal 
and  flexural  modes*.  Thus  the  bars  could  be  exoited  in  simple  axial  vibration 
and  their  natural  frequencies  measured  as  a function  of  the  external  pressure* 
Results  of  such  experiments  are  ah  ova  in  Pig.  12.  Two  sets  of  data  are  given  t 
one  with  bells  having  a dimensional  tolerance  of  t 50  x 10 in.  and  the 
other  1 10  x 10~^  In*  As  may  be  seen,  the  frequencies  of  the  bar  made  with 
the  better  balls  are  oloeer  to  the  theoretical  frequencies  and  the  agreement 
Improves  In  both  cases  with  increasing  pressure*  The  reason  for  this  becomes 
apparent  when  the  dimensional  tolerances  are  compared  with  the  relative 
approach  of  the  balls  under  the  initial  pressure.  When  <f»  =*  2 psi, 

* a 1.955  x 10"6  in.  and  when  a 14*7  psi,  °t  3 7.39  x 10“^  in.  Thus 
many  spheres  may  be  expected  to  be  under  larger  and  smaller  initial  contact 
forces  than  if  all  spheres  were  identical  in  size  and,  also,  some  spheres 

s 

may  be  loose*  It  may  be  shown  that  the  presence  of  off -size  or  loose  spheres 
diminishes  the  stiffness  (and  hence  the  frequency  of  vibration)  of  the  array 
and  the  diminution  becomes  greater  with  increased  spread  of  the  dimensional 
tolerance  and  reduction  of  pressur-. » These  effects  are  reflected  in  the  data 
shown  in  Fig.  12* 

Measurements  of  logarithmic  decrement  of  the  vibrations  were  also  made, 
but  they  cannot  be  compared  with  the  theory  until  the  imaginary  parts  of  the 
compliances  are  introduced  into  Equations  (14)* 
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Fig,  8s  Energy  dissipation  per  cycle  as  a 
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with  experimental  data  by  Johnson, 
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